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ABSTRACT 


The  shielding  requirements  for  the  protection  of  the  crew  of  a  manned 
orbital  command  post  against  the  natural  space  radiations  are  investigated. 

Two  types  of  orbits  of  military  importance  and  of  wide  applicability  are 
considered — a  long-duration,  high-altitude  orbit  above  the  Van  Allen  Belt  and 
a  short-duration,  low-altitude  polar  orbit  below  it.  Model  environments  for 
both  orbits  in  terms  of  solar  flare,  cosmic,  and  Van  Allen  Belt  radiations  are 
postulated.  Radiobiological  tolerance  criteria  are  investigated,  and  a  some¬ 
what  unique  criteria,  based  on  partial  recovery  of  sustained  somatic  damage, 
is  proposed  for  the  long  duration  mission. 

A  mathematical  description  of  the  radiation  transport  of  the  separate 
environmental  components  through  the  radiation  shield  is  formulated.  Appro¬ 
priate  simplifications  are  used  to  obtain  expressions  for  the  doses  due  to 
primary  protons,  secondary  protons  and  neutrons,  and  bremsstrahlung.  Calcu¬ 
lations  are  carried  out  on  the  IBM  7030  computer  to  obtain  dose  vs.  thickness 
curves  for  different  types  and  amounts  of  shielding  material.  These  curves, 
together  with  the  assumed  model  environment,  and  the  postulated  radiobiological 
tolerance  criteria,  are  used  to  calculate  minimum  shielding  thicknesses  for 
both  types  of  orbits.  Conclusions  are  then  drawn  on  the  total  amount  of  radiation 
shielding  material  that  must  be  carried  for  each  orbit. 
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SECTION  I 


INTRODUCTION 

It  is  vital  to  our  national  defense  posture  to  always  maintain  effective 
command  and  control  over  our  global  strike  forces.  This  requires  that  we 
have  command  posts  that  are  survivable  in  the  event  of  a  surprise  attack.  As 
pointed  out  by  General  Power,  a  promising  method  of  attaining  such  surviv¬ 
ability  is  through  the  use  of  maneuverable  command  posts  in  space. 

Over  the  past  several  years,  the  merits  of  various  t3'pes  of  manned  orbital 

command  posts  (MOCOPs)  have  been  investigated.  The  vulnerability  of  such 

[2] 

vehicles  to  direct  nuclear  attack  has  been  reported  previously.  In  investigat¬ 
ing  the  feasibility  of  such  vehicles,  one  of  the  most  important  considerations  is 
the  protection  required  to  make  them  survivable  in  the  natural  space  environ¬ 
ment. 

Perhaps  the  most  important  component  of  the  space  environment  for 
manned  space  flight  is  the  charged  particle  radiations.  These  radiations  may 
be  divided  into  the  three  following  t3q)es:  Van  Allen  Belt  radiations  (fjrotons 
and  electrons  trapped  in  the  geomagnetic  field);  solar  flare  radiations  (transient, 
energetic  proton  streams,  associated  with  flares  on  the  sun);  and  cosmic  radia¬ 
tions  (protons  and  heavier  ions  which  arrive  from  all  directions  of  our  galaxy). 
Because  of  the  biological  effects  produced  on  the  crew  by  tliese  radiations, 
shielding  must  be  employed  in  the  space  vehicle  to  attenuate  the  radiations.  For 
cer  tain  types  of  orbits  and  boosters  available  in  a  given  time  period,  the  amount 
ot  shielding  required  to  give  even  minimum  protection  to  the  crew  may  so  severely 
limit  the  payload  as  to  render  the  mission  infeasible.  For  tliis  reason, 


1 


evaluation  of  the  merit  and  feasibility  of  proposed  types  of  command  and  control 
missions  in  space  necessitates  careful  investigation  of  radiation  protection 
requirements. 

The  general  problem  of  the  vulnerability  of  MOCOPs  to  space  radiations 
is  investigated  by  considering  two  specific  tjqjes  of  orbits.  The  first  is  a  high- 
altitude  (50,  000  n.  m. )  orbit  and  the  second  is  a  low-altitude  (200  n.  m. )  polar 
orbit.  These  orbits  were  chosen  both  because  of  their  military  significance 
and  because  each  one  is  representative  of  a  class  of  radiation  vulnerability 
problem.  The  results  that  are  obtained  in  this  study  will  be  applicable  to  a 
fairly  broad  range  of  orbits.  For  instance,  many  of  the  results  for  the  high- 
altitude  orbit  will  be  applicable  to  any  orbit  in  cislunar  space  that  is  outside  the 
Van  Allen  Belt  most  of  the  time,  and  even  to  interplanetary  transfers.  The 
results  for  the  low-altitude  polar  orbit  will  be  applicable  to  any  polar  orbit 
below  the  Van  Allen  Belt  and  will  be  an  upper  bound  on  any  nonpolar  orbit 
below  the  Van  Allen  Belt. 

This  paper  first  discusses  the  assumed  mission  profiles  and  vehicle 
configuration.  The  environment  models  assumed  for  the  missions  are  then 
presented.  Radiation  dose  criteria  for  the  crew  are  next  stipulated.  'ITie 
shielding  calculations  are  briefly  discussed  and  are  presented  in  detail  in 
Appendices  I  and  11.  Finally,  the  results  of  the  shielding  calculations  and  the 
dose  criteria  are  used  to  determine  requisite  amounts  of  shielding. 
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SECTION  n 


MISSION  PROFILES  AND  VEHICLE  CONFIGURATION 

HIGH-ALTITUDE  ORBIT  MISSION  PROFILE 

A  high-altitude  orbit  for  the  MOCOP  is  attractive  because  it  is  difficult 
for  an  enemy  to  successfully  launch  a  surprise  attack  on  it,  and  also  because 
the  MOCOP  has  considerable  maneuvering  room  to  evade  enemy  detection  and/or 
attack.  The  orbit  should  not  be  so  high  as  to  make  contact  with  the  earth 
difficult  but  should  be  high  enough  to  be  outside  of  the  trapped  radiation  belt. 

With  these  considerations  in  mind,  an  operational  altitude  of  50,  000  nautical 
miles  was  selected.  As  the  results  of  the  study  later  showed,  the  choice  of 
particular  orbital  altitude  is  not  too  important  in  determining  radiation  shield¬ 
ing  requirements,  as  long  as  it  is  outside  the  geomagnetic  field.  The  radiation 
encountered  in  this  orbit  is,  as  far  as  is  known,  the  same  as  that  encountered 
everywhere  in  cislunar  space  and  will  be  designated  as  the  "free-space" 
radiation. 

It  is  assumed  that  the  vehicle  arrives  at  the  operational  altitude  bv  first 
going  into  a  100  u.  m.  parking  orbit  and  then  transferring  to  the  higher  altitude 
by  means  of  a  Hohmann  semi-ellipse.  Once  the  vehicle  reaches  operational 
altitude,  it  is  assumed  that  it  is  continuously  manned  by  rotating  crews.  For 
this  type  of  orbit,  it  is  desirable  to  keep  crew  rotations  to  a  minimum.  There¬ 
fore,  extended  orbital  tours  of  duty  in  the  range  of  three  months  to  one  year 
will  be  considered. 
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LOW-ALTITUDE  ORBIT  MISSION  PROFILE 


While  a  low-altitude  orbit  for  a  MOCOP  may  be  more  vulnerable  to  enemy 
attack,  it  does  possess  several  attractive  features.  For  instance,  such  orbits 
have  less  demanding  initial  and  resupply  booster  requirements  than  high-altitude 
orbits,  allow  observation  of  the  earth  in  greater  detail,  and,  as  will  be  shown, 
can  permit  less  stringent  radiation  shielding  requirements.  (It  is  assumed  that 
we  are  considering  a  time  period  in  which  the  geomagnetic  field  at  low  altitudes 
is  essentially  free  of  Argus-type  electrons. )  For  MOCOP  missions,  polar  and 
nearly-polar  orbits  are  of  interest  because  of  their  earth-coverage  properties. 
Hence,  the  low-altitude  orbit  that  is  considered  is  the  polar  circular  orbit  at 
a  nominal  altitude  of  200  n.  m. 

Because  of  the  relative  ease  of  resupply  and  crew  rotation  as  compared 
to  the  high-altitude  case,  it  is  envisioned  that  the  orbital  tours  of  duty  will  be 
shorter.  Therefore,  tours  of  duty  in  the  range  from  one  week  to  three  months 
will  be  considered. 

VEHICLE  CONFIGURATION 

For  both  the  high-  and  low-altitude  cases,  the  MOCOPs  are  envisioned  as 
one  or  several  identical,  semi-independent,  rigid  modules.  Each  module  is 
assumed  to  be  cylindrical  in  shape  with  a  diameter  of  12  feet  and  a  length  of  25 
feet.  The  effects  of  the  radiations  in  a  typical  module  are  considered  by  neglect¬ 
ing  the  interaction  between  the  modules,  which  is  estimated  to  be  small.  It  is 
further  assumed  that  the  modules  have  no  reentry  capability  in  themselves. 

Each  module  consists  of  an  outer  cylindrical  shell  plus  an  inner,  spherical 
"storm  cellar"  for  protection  against  severe  radiations. 
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The  design  of  the  outer  shell  is  influenced  by  launch  and  orbit  structural 

loads  as  well  as  micrometeorite  protection  and  thermal  balance  considerations. 

It  is  assumed  that  the  shell  is  of  aluminum-foamed  polyurethane-aluminum 

2  2 

multiwall  construction  with  a  unit  weight  of  1.  5  gm/cm  («  3.  1  Ib/ft  ).  The 
storm  cellar  concept  for  radiation  shielding  is  used  to  minimize  the  weight  of 
shielding  material  and  the  storm  cellar  is  occupied  only  in  times  of  severe  solar 
flares  and,  in  the  high-altitude  orbit  case,  passage  through  the  Van  Allen  Belt. 
Since  the  time  between  the  observance  of  the  maximum  visual  flare  intensity  and 
the  arrival  of  substantial  proton  fluxes  near  the  earth  is  greater  than  a  half- 
hour,  there  is  adequate  time  for  the  crew  to  retire  to  the  storm  cellar. 


SECTION  m 


MODEL  ENVIRONMENT 

HIGH-ALTITUDE  ORBIT  ENVIRONMENT 
Solar  Flares 

The  current  state  of  knowledge  concerning  the  solar  flare  phenomenon 
makes  the  postulation  of  a  model  solar  flare  environment  difficult  because  of 
the  following  basic  deficiencies  in  our  knowledge.  First,  there  is  little  under¬ 
standing  of  the  mechanism  of  the  genesis  of  a  solar  flare  and  its  relation  to 
associated  proton  streams.  Their  occurrence  seems  to  follow  a  random  pattern 
and  each  flare  appears  to  have  its  own  characteristic  time-varying  energy 
spectrum.  Second  the  limited  number  of  observations  rule  out  any  meaningful 
statistical  analysis  of  flares.  It  has  long  been  known  that  solar  activity  follows 
an  11 -year  cycle  and  recent  observations  indicate  that  most  important  solar  flares 
occur  before  and  after  the  maximum  of  the  cycle.  However,  the  most  important 
data  on  flares  have  only  been  obtained  through  approximately  one  solar  cycle. 
Third,  even  fairly  recent  observations  are  often  lacking  in  data  at  low  energies 
and  early  times.  Hopefully,  in  the  future  these  three  deficiencies  will  be 
rectified  by  continued  observations  and  utilization  of  data  from  such  special- 
purpose  satellites  as  OSO  and  IMP.  For  the  present,  one  is  forced  to  rely  on 
past  observations  to  hypothesize  what  solar  flare  environment  will  be  encountered 
by  future  astronauts. 

The  important  solar  flares  that  can  deliver  a  lethal  or  disabling  sublethal 
dose  to  a  space  crew  are  sometimes  categorized  as  "major"  and  "giant"  events. 

A  major  solar  flare  has  a  high-flux,  low-to-medium  energy  spectrum.  Typical 
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of  the  most  intense  of  this  type  are  the  May  10,  1959  flare  (high  flux,  low 

energy)  and  the  November  12,  1960  flare  (high  flux,  medium  energy).  Major 

flares  occur  two  to  four  times  a  year  during  periods  of  increasing  or  decreasing 

solar  activity;  e.  g. ,  1959  was  possibly  the  worst  year  for  major  solar  flares 

as  four  occurred  that  would  have  caused  an  integrated  skin  dose  in  free  space  of 

2  [3  4l 

greater  than  100  rad  behind  1  gm/cm  of  shielding.  ’  Giant  solar  flares 
have  low  or  moderate  fluxes  of  extremely  high  energy  (in  the  relativistic  range), 
and  occur  maybe  once  every  two  years.  The  classic  example  of  this  type  is  the 
February  23,  1956  flare,  the  most  energetic  ever  observed. 

Other  flares  of  lesser  flux  and  energy  occur  much  more  frequently  and  it 
is  conceivable  that  their  cumulative  effect  may  be  important.  However,  on  the 
basis  of  estimates  of  Reference  3,  even  if  the  crew  remains  outside  of  the 
radiation  shield  during  minor  flares,  they  will  accumulate  a  dose  of  no  greater 
than  12  rad  in  a  year  from  these  flares.  For  the  present  study,  it  is  assumed 
that  the  crew  uses  the  radiation  shield  as  a  shelter  against  the  more  important 
of  these  minor  flares  so  that  the  cumulative  dose  is  reduced  to  a  very  low  value 
which  can  be  neglected. 

The  solar  flare  model  environment  postulated  must  be  somewhat  arbitrary, 
and  hence  given  to  second-guessing.  In  the  last  analysis,  it  must  be  arrived  at 
by  judgment,  balancing  conservatism  against  probability.  There  is  no  apparent 
physical  reason  why  six  giant  flares  such  as  that  of  February  23,  1956  cannot 
occur  during  a  twelve  month  period,  but  on  the  basis  of  past  experience  such  an 
occurrence  seems  very  unlikely.  The  present  study  considers  two  types  of 
flares:  Type  A,  a  high-flux,  low-energy  flare  that  is  represented  by  the  May  10, 
1959  flare;  and  Type  B,  (in  reality  two  consecutive,  high-flux,  medium-energy 
flares)  that  is  represented  by  the  November  12-15,  1960  sequence.  This  latter, 
double-type  flare  is  considered  because  the  radiation  damage  criteria  to  be  used 
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is  more  sensitive  to  two  identical  flares  occurring  consecutively  than  to  two 
identical  flares  occurring  several  months  apart.  These  particular  flares  were 
chosen  for  the  model  environment  because  they  are  among  the  most  intense  of 
their  types  and  because  fairly  good  data  exists  on  their  behavior. 

Orbital  tours  of  duty  of  from  three  months  to  one  year  are  considered  and 
the  assumed  time-phasing  of  the  flares  is  shown  below  for  3-,  6-,  and  12-month 
missions. 

12  3  456  789  10  11  12 

12  month  mission  ABA 

6  month  mission  A  B 

3  month  mission  B 

A  =  May  10,  1959  flare 
B  =  November  12-15,  1960  flare  sequence 

The  environment  is  specified  in  terms  of  A-  and  B-type  flares,  but  the 
results  will  be  analyzed  to  obtain  the  effect  of  substituting  a  February  23,  1956 
type  flare  for  an  A  or  B  type. 

Having  selected  the  flares  to  be  used  in  the  model  environment,  the 

problem  arises  as  how  best  to  represent  the  flare  behavior  analytically.  After 

some  consideration,  it  was  decided  to  use,  wherever  possible,  the  exponential 

[5  3] 

rigidity  formulation  proposed  by  Freier  and  Webber.''  ’  This  is  a  one 
parameter  spectrum  of  the  form 


-R/R  (t) 

J  =  J  (t)  e 
o 


(1) 
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where 


J  =  integral  proton  flux  with  rigidity  greater  than  R 
(protons/cm^  sec  ster) 

R  =  proton  rigidity  (momentum  per  unit  charge) 

=  total  integral  proton  flux  (R  >  0) 

R  =  characteristic  rigidity  (for  a  given  time  in  a  given  event) 
o 

Specific  analytical  representation  of  the  May  10,  1959,  November  12-15, 

1960,  and  the  February  23,  1956  flares  are  discussed  below. 

November  12-15,  1960  Flare  Sequence 

Values  of  (t)  and  R^  (t)  are  taken  directly  from  Reference  5  and  are 
presented  in  Table  1. 

May  10,  1959  Flare 

Values  of  J  (t)  and  R  (t)  from  t  +32  hours  onward,  where  t 

o  o  o.  m.  o.  m. 

is  the  time  of  optical  maximum,  are  given  in  Table  1'  .  Values  of  R^  (t) 

between  t  =  t  and  t  =  t  +32  hours  are  obtained  by  extrapolation 

o.  m.  o.  m. 

of  the  tabulated  values  of  R  for  t  >  32  hours.  Values  of  J  (t)  between 

o  o 

t  =  t  and  t  =  t  +32  hours  are  obtained  by  assuming  that  J  rose 

o.m.  o.m.  o 

linearly  for  20  hours  (based  on  data  in  References  3  and  4)  until  it  reached  the 

value  of  J  at  t  =  t  +32  hours.  J  is  taken  constant  between  20  and 
o  o.m.  o 

32  hours. 

February  23,  1956  Flare 

Analytic  representation  of  this  flare  is  hampered  by  lack  of  observational  data. 
Sufficient  data  does  not  exist  to  represent  this  flare  in  terms  of  exponential  rigidity 
spectra  so  that  a  formulation  commensurate  with  available  data  is  necessary. 
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Table  1 


J  and  R  Values  for  November  12-15,  1960  and  May  10,  1959  Flares* 
o  o 


Date  of  Solar  Flare 

Time  of 
Optical 
Maximum 

Time  of 

Measurement 

J  ,  protons/ 

0 

2 

c  sec  ster 

m 

R  ,  Mv 
o 

May  10,  1959 

2118 

5/12  0500 

30  X  10^ 

65 

0520-, 

0915^ 

3 

18  X  10 

70 

0915-, 

1335^ 

4. 1  X  10^ 

70 

1335-. 

1750^ 

2. 8  X  10^ 

66 

November  12,  1960 

1329 

11/12  1930 

1100 

280 

2000 

1250 

240 

11/13  0200 

3400 

185 

0800 

3400 

155 

1305 

3000 

120 

1830-. 

2000^ 

2800 

95 

2230 

1000 

105 

11/14  0500 

1000 

95 

November  15,  1960 

0221 

11/15  0500 

250 

375 

1130 

320 

240 

1030-, 

1230^ 

285 

175 

2130 

3300 

120 

11/16  0930 

1000 

100 

1430-, 

1730^ 

300 

100 

*  (From  Table  2  of  Reference  5) 
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[3  4] 

It  is  assumea  ’  ^  that  for  a  given  energy  the  timewise  variation  in 
integral  proton  flux  is  given  by 


J  (t,  E  >  E  \  =  J  (E  >  EJ  e 
'  1'  max  '  1' 


-t\  (Ej) 


(2) 


where 

J  (E  >  E  1  =  timewise  maximum  integral  flux 

max  '1' 

t  =  time  measured  from  maximum  flux 

tj^  (E^)  =  characteristic  decay  time  (which  is  a  function  of 

energy) 

Integral  energy  flux  spectra  (J  versus  E)  are  given  by  Foelsche^^^  for 
prompt  and  19-hour  conditions.  Because  of  the  extreme  energy  of  this  flare, 
the  rise  time  is  short,  so  the  assumption  is  made  that  the  rise  time  is  zero, 
and  that  the  prompt  spectrum  represents  the  maximum  flux  condition.  Then, 
with  the  prompt  and  19-hour  spectra  data  from  Reference  6,  Eq.  (2)  can  be  used 
to  obtain  t^^  (E). 

It  is  assumed  that  the  differential  energy  spectrum  follows  a  piecewise  power 
th 

law,  and  for  the  i  interval  has  the  form 


-t/t^  (E) 


D 


-n. 


(E) 


C  e 
1 


(E.  <  E  <  E  ) 
'  i  i+1^ 


(3) 


where 

=  differential  proton  flux  in  i*'^  interval 
th 

C.,  n.  =  constants  in  i  interval 


The  integral  proton  flux  can  be  obtained  by  integration  of  Eq.  (3).  At  t  =  0 
this  is 


J  (E  >  E.)  =  J  (E  >  E.^^)  + 

One  can  then  also  write 


J  (E  >  E.)  -  J  (E  >  E^^^) 
J  (E  >  E.)  -  J  (E  >  E.^^) 


where  E.  <  E.  <  E.  , 

1  J  1+1 

Values  of  n^  can  be  obtained  from  a  given  integral  energy  prompt  spectrum  by 

reading  off  values  of  J  and  E  at  the  end  points  (corresponding  to  E^  and  1-.^^) 

and  an  intermediate  point  (corresponding  to  E^)  of  the  various  intervals,  and 

then  applying  Eq.  (5).  With  m  determined,  Eq.  (4)  may  be  used  to  calculate 

the  various  C.. 

1 

The  procedure  just  indicated  is  applied  to  the  data  of  Reference  G  to  calcu¬ 
late  the  quantities  C^,  t^^  ,  and  n^  in  Elq.  (3)  that  define  the  spectrum.  The 

results  are  given  in  Table  2.  It  should  be  emphasized  that  because  of  the 
scarcity  of  data  on  this  flare,  and  the  resultant  assumptions  and  extrapolations 
that  were  necessary  to  represent  its  behavior,  the  confidence  level  in  the 
representation  is  lower  than  for  the  A-  and  B-type  flares. 
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Table  2 


Constants  for  February  23,  1956  Flare 


Interval 

Number 

Energy  Range,  Mev. 

^D. ,  hr 

1 

n 

1 

_2  -1  -1  ”i"^ 

C  ,  cm  sec  ster  Mev 

1 

1 

20  - 

100 

4.  8 

1.  0 

4.  66  X  10^ 

2 

100  - 

160 

4.4 

1.  0 

4.  68  X  10^ 

3 

160  - 

300 

3.  6 

1.  08 

5.  14  X  10^ 

4 

300  - 

550 

2.  9 

2.  16 

6 

2.  75  X  10 

5 

550  - 

1,  600 

2.  3 

2.  37 

3.  06  X  10^ 

6 

1, 600  - 

5,  000 

1.  65 

5.  21 

7.  57  X  10^® 

7 

5, 000  - 

10, 000 

1.  65 

7.  59 

3.  93  X  10^^ 

Van  Allen  Belt 

It  is  assumed  that  the  traversal  of  the  Van  Allen  Belt  from  parking  in  high- 
altitude  orbit,  and  vice  versa,  is  made  in  the  plane  of  the  geomagnetic  equator. 
Orientation  of  the  transfer  trajectory  in  this  plane  will  result  in  a  maximum 
dose  and  thus  furnish  an  upper  limit  on  the  trapped  particle  dose.  A  Hohmann 
semi-ellipse  is  chosen  for  the  transfer  trajectory,  giving  a  transit  time  of  16.  7 
hours.  The  result  of  calculations  carried  out  to  determine  the  vehicle's  altitude 
as  a  function  of  time  during  the  transfer  is  shown  in  Fig.  1. 

Proton  Flux 

Several  proposed  curves  of  differential  proton  energy  spectra  in  the  heart 

of  the  inner  belt  are  given  in  Fig.  3-2  of  Reference  7.  The  recommendation 

(8) 

therein  to  use  the  curve  based  on  Lockheed  HEPDEX  and  LEPDEX  data  is 
followed:  This  spectrum  is  defined  as 
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Fig.  1 .  Altitude  vs  Time  During  Transfer  From  Parking  to  High-Altitude  Orbit 


3  2 

(p  =  5.  33  X  10  protons/cm  -sec-Mev  E  <  30  Mev 

(6) 

8  —3  0  2 

(p  =  1.44  X  10  E  ’  protons/cm  -sec-Mev  E  >  30  Mev 

For  the  present  investigation,  it  is  assumed  that  this  spectrum  is  invariant 

with  altitude.  A  curve  of  integrated  proton  flux  (E  >  30  Mev)  versus  altitude, 

[8] 

based  on  Explorer  IV,  Pioneer  III,  and  HEPDEX  data  from  Fig.  3-3  of 
Reference  7  is  shown  in  Fig.  2. 

Note  that  the  dose  rate  behind  a  shield  is  proportional  to  (p,  and  conse¬ 
quently  to  J.  Let  the  dose  rate  behind  a  given  shield  in  the  heart  of  the  inner  belt 
be  denoted  by  D  ).  Then  at  any  other  altitude: 
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PftOTON  rtux.  E  >  30  liEV,  mOTDNS/CM*  sec 


Fig.  2.  Integrated  Proton  Flux  in  the 
Van  Allen  Belt  vs  Altitude 


D(h) 


J  (h,  E  >  30) 

J  E  >  30) 


(7) 


where 


J  (h,  E  >  30) 

J  E  >  30) 


reference  altitude  =  altitude  at  peak  of  inner  belt 

integral  flux  at  altitude  h  as  obtained  from  Fig.  2 

integral  flux  at  reference  altitude  as  obtained 
from  Fig.  2. 


The  total  proton  dose  D  received  during  transit  is 


D 


D  (t)  dt  =  D  (h^^^) 


J  (h,  E  >  30) 

J  (h  E  >  30) 
ref 


(8) 


16 


where 


T  =  transit  time 

One  can  combine  the  curve  of  h  versus  t  from  Fig.  1  with  the  curve  of 
J  (E  >  30)  versus  h  in  Fig.  2,  to  obtain  the  curve  of  J  (E  >  (E  >  30) 

versus  t  in  Fig.  3.  Numerical  integration  of  this  curve,  in  accordance  with 
Eq.  (8),  then  leads  to 

D  =  T.  D  (9) 

where  T*  is  the  integral  in  the  right-hand  expression  in  Eq.  (8).  Equation  (9) 

may  be  interpreted  as  the  proton  dose  received  during  the  transfer  trajectory 

equals  the  product  of  the  dose  rate  received  in  the  heart  of  the  inner  belt  and  a 

fictitious  time  T*. 

P 


Fig.  3.  Nondimensional  Integrated 
Proton  Flux  in  the  Van 
Allen  Belt  vs  Transfer  Time 
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The  above  formulation  of  the  problem  requires  only  one  shielding  calcula¬ 
tion  for  a  given  shield  configuration.  This  calculation  is  carried  out  for  the  peak 
flux  and  the  result  is  multiplied  by  the  factor  T*  =  0. 128  hour  to  yield  the 
proton  dose  during  orbit  transfer.  For  other  transfer  trajectories,  the  same 
computational  procedure  can  be  followed  with  appropriate  values  of  T*. 

Electrons 

The  undisturbed  electron  spectra  encountered  by  the  vehicle  in  passing 
through  the  Van  Allen  Belt  is  formulated  in  a  manner  similar  to  that  used  for 
the  Van  Allen  proton  spectra  previously  described.  In  this  case,  it  is  assumed 
that  there  are  two  invariant  spectra — one  for  the  inner  belt  and  one  for  the  outer 
belt.  The  assumed  spectra  are  shown  in  Fig.  4,  as  taken  from  Reference  9,  the 


Fig.  4.  Integral  Electron  Flux  Spectra  in  the  Van  Allen  Belt 
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inner  belt  spectra  being  that  given  in  Reference  10  and  the  outer  belt  that  given 
by  Van  Allen  in  the  symposium,  Space  Flight  Report  to  the  Nation.  Also  presented 
in  Reference  9  are  data  from  which  it  is  possible  to  construct  the  curve  of 
integral  electron  flux  (E  >  20  kev)  versus  altitude  that  is  shown  in  Fig.  5. 

Preliminary  shielding  calculations  showed  that  the  dose  due  to  the  elec¬ 
trons  could  be  ignored  compared  to  that  due  to  their  bremsstrahlung,  so  that 
only  the  latter  is  considered  here. 

Calculation  of  the  bremsstrahlung  dose  is  carried  out  in  a  manner  analogous 
to  that  described  in  the  previous  section  for  the  protons,  except  that  separate 
integrations  are  carried  out  over  the  inner  and  outer  belts.  Therefore,  analogous 
to  Eqs.  (8)  and  (9),  there  are 


(10) 


(11) 


where  the  subscript  e  denotes  electron,  and  the  subscripts  i  and  o  denote 

inner  and  outer  belt,  respectively.  Figures  4  and  5  can  then  be  combined  to 

obtain  the  curve  of  J  /J  shown  in  Fig.  6.  Numerical  integration  of  Fig.  6 

ref 

yields  T*  .  =  0.919  hour  and  T*  =  0. 607  hour, 
e,  1  e,  o 
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(E>20KEV  ) 


Fig.  5.  Integrated  Electron  Flux  in 
the  Van  Allen  Belt  vs 
Altitude 


Fig.  6.  Nondimensional  Integrated 
Electron  Flux  in  the  Van 
Allen  Belt  vs  Transfer  Time 
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Cosmic  Rays 


It  is  assumed  that  except  during  periods  of  major  and  giant  solar  flares, 
the  normal  on-orbit  operating  mode  of  the  crew  is  outside  the  storm  cellar. 
Hence,  during  most  of  their  duty  tour,  the  crew  will  be  subjected  to  the  effects 
of  the  background  cosmic  radiation  (or  galactic  cosmic  radiation)  shielded  only 
by  the  wall  of  the  space  vehicle.  These  radiations  are  of  such  high  energy  that 
the  vehicle's  wall  will  not  significantly  attenuate  them,  and  the  dose  received 
inside  the  cabin  will  be  essentially  unshielded. 

While  the  spectrum  and  composition  of  the  cosmic  rays  are  thought  to  be 

known  fairly  well  (e.g.  ,  References  11,  12),  the  biological  effects  of  the  high- 

energy  protons  and  the  heavy  primaries  are  at  present  poorly  understood.  In 

view  of  this  latter  uncertainty,  it  is  inconsistent  to  use  a  detailed  cosmic  ray 

spectrum  in  the  shielding  calculations.  Instead,  the  unshielded  dose  rate  of 

[13] 

26  rem/yr  proposed  by  Tobias^’  ^  will  be  accepted. 

LOW-ALTITUDE  ORBIT  ENVIRONMENT 
Geomagnetic  Field  Effects 

The  only  radiation  components  encountered  by  a  vehicle  in  the  low-altitude 
orbit  are  those  of  solar  flares  and  cosmic  rays,  since  the  orbital  altitude  is 
below  the  Van  Allen  Belt.  The  main  difference  in  radiation  environment  encoun¬ 
tered  in  this  orbit  and  that  encountered  in  the  high-altitude  orbit  is  due  to  the 
filtering  action  of  the  geomagnetic  field  on  charged  particles.  *  This  effect 


■•There  is  another  difference  between  the  high-  and  low-altitude  radiation 
environments  that  was  not  considered  in  the  present  study  because  of  lack  of 
information,  and  because  its  neglect  errs  on  the  conservative  side.  This  differ¬ 
ence  is  that,  due  to  the  shielding  action  of  the  earth,  the  incident  radiations  may 
not  be  isotropic  in  low-altitude  orbits.  It  is  not  possible  to  estimate  the  import¬ 
ance  of  this  effect,  but  an  upper  bound  is  that  it  reduces  the  flux  by  one  half. 
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is  a  function  of  geomagnetic  latitude  and  for  each  latitude  there  is  a  correspond¬ 
ing  minimum  rigidity  (cutoff  rigidity)  that  a  charged  particle  must  possess  to 
penetrate  this  field.  The  cutoff  rigidity  varies  from  about  14.9  Bev/C  (see 
Appendix  II)  at  the  geomagnetic  equator  to  what  is  thought  to  be  about  zero  at 
the  geomagnetic  poles  for  normally  incident  particles.  Thus,  a  vehicle  in  an 
orbit  in  the  plane  of  the  geomagnetic  poles  receives  a  variable  dose  rate,  even 
if  the  free-space  dose  rate  is  constant.  The  dose  rate  during  orbit  will  be  the 
free-space  dose  rate  when  the  vehicle  is  over  either  of  the  geomagnetic  poles 
and  close  to  zero  when  it  is  in  the  equatorial  regions.  For  a  given  (low)  altitude, 
such  a  vehicle  will  receive  a  maximum  dose  rate  as  compared  to  vetiicles  in 
orbits  with  different  inclinations.  Since  the  geomagnetic  and  geographic  poles 
are  fairly  close  (the  geomagnetic  north  pole  is  at  78.  5^  N  latitude),  the  dose 
received  during  a  geomagnetic  polar  orbit  should  furnish  a  fairly  close  and 
conservative  estimate  of  that  received  during  a  geographical  polar  orbit. 

Because  of  the  much  greater  computational  simplicity  in  calculating  the  radia¬ 
tion  dose  in  a  geomagnetic  polar  orbit,  as  compared  to  geographical  polar 
orbit,  the  radiation  environment  assumed  for  the  low-altitude  orbit  is  taken  to 
be  that  in  the  plane  of  the  geomagnetic  poles. 

Solar  Flares 

The  solar  flares  that  are  considered  are  the  same  ones  considered  for  the 
high-altitude  orbit,  i.e.,  May  10,  1959,  November  12-15,  1960,  and  February 
23,  1956.  The  free-space  spectra  of  those  flares  are  given  at  the  beginning  of 
Section  3.  The  manner  in  which  the  earth's  magnetic  field  is  taken  to  modify 
the  free-space  spectra  is  discussed  in  Appendix  II. 
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For  the  model  environment,  it  is  assumed  that  during  the  mission,  there 
is  an  occurrence  of  one  of  the  three  types  of  flares.  Because  a  cumulative 
radiobiological  tolerance  criterion  is  assumed,  the  time  during  the  mission 
when  the  flare  occurs  is  immaterial. 

Cosmic  Rays 

In  order  to  evaluate  the  effects  of  the  geomagnetic  field  on  the  dose  due 
to  the  cosmic  rays,  it  is  necessary  to  look  at  the  cosmic  ray  spectrum  in  more 
detail  than  was  done  for  the  high-altitude  orbit.  A  cosmic  ray  spectrum  of  the 
following  kind  is  assumed: 

for  4  <  E  <  380  Mev: 

(E)  =  0.  514  X  lO”^  (1  +  lO”^  E)®’ 

for  380  ^  E  <  oo  Mev:  (12 1 

N  (E)  =  19.  1  X  10“^  (1  +  lO"®  E)“^’  ^ 

4b 


where 

2 

N  and  N  =  differential  proton  flux  in  particles/(cm  -sec-ster-Mevl 

-L  ^ 

E  =  proton  energy  in  Mev 

The  spectrum  given  by  Eq.  (12)  was  based  on  data  given  in  Fig.  30  of  Reference 
12,  and  normalized  to  agree  with  the  free-space  dose  of  26  rem/yr  that  was 
assumed  above. 

The  effect  of  the  geomagnetic  field  on  the  cosmic-ray  dose  is  essentially 
the  same  as  on  the  solar-flare  dose,  except  that  the  free-space  cosmic-ray 
dose  is  time-independent.  The  calculation  of  tlie  cosmic-ray  dose  for  the  low- 
altitude  orbit  is  discussed  in  Appendix  II. 
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SECTION  IV 


RADIOBIOLOGICAL  TOLERANCE  CRITERIA 

GENERAL  CONSIDERATIONS 

The  ideal  approach  to  the  radiation  shielding  problem  is  to  stipulate  the 
worst  conceivable  radiation  environment  and  design  a  shield  behind  which  the 
crew  experiences  no  discernible  radiobiological  effects.  Such  a  shield  would 
be  very  massive,  and,  unfortunately,  the  limited  capabilities  of  boosters  of  the 
present  and  foreseeable  future  precludes  such  an  approach.  Thus  weight 
considerations  force  us  to  be  a  little  less  conservative  in  stipulating  the  environ¬ 
ment,  and  a  little  more  tolerant  of  minor  biological  effects.  In  the  previous 
sections,  we  postulated  a  radiation  environment  that  appears  reasonable  based  on 
present  data  and  slightly,  but  not  overly,  conservative.  We  will  now  consider 
what  constitutes  tolerable  radiobiological  effects. 

The  military  has  traditionally  been  asked  to  face  higher  risks,  in  perform¬ 
ance  of  their  duties  involving  national  security,  than  have  civilians.  Extending 
this  reasoning  a  little  further,  personnel  on  military  space  operations  could 
conceivably  be  asked  to  face  slightly  greater  risks  than  personnel  on  civilian 
scientific  or  exploratory  space  missions.  In  the  context  of  military  space 
operations,  a  risk  might  also  involve  the  possibility  of  some  biological  effects  of 
space  radiations  that  manifest  themselves  long  after  the  flight  is  over,  e.  g.  , 
general  life-shortening.  The  maximum  risk  that  the  crew  should  be  asked  to 
face  should  be  a  function  of  the  importance  of  the  mission  to  national  security 
(e.  g. ,  a  routine  cold-war  mission  should  entail  less  risk  than  a  general-war 
mission).  Also,  in  discussing  radiation  criteria,  one  should  consider  not  only  long¬ 
term  or  chronic  effects  but  also  the  acute  effects  that  occur  during  the  flight. 
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For  example,  while  obviously  the  crew  should  not  sustain  a  lethal  dose,  they  also 
should  not  sustain  a  temporarily  disabling  dose  that  will  decrease  their  efficiency 
below  that  necessary  to  j)erform  their  mission.  The  magnitude  of  this  latter 
tyi)c  of  dose  is  a  function  of  the  mental  and  physical  difficulty  of  the  task  to  be 
performed.  Fxecution  of  command  is  a  relatively  sophisticated  task  and  its 
performance  will  not  tolerate  much  degradation  in  the  commander's  normal 
efficiency,  thus  implying  that  the  commander  should  only  be  subjected  to  a 
relatively  low  dose.  Another  facet  of  the  problem  that  should  be  considered  is 
that  the  effect  of  a  large  dose  received  during  the  mission  might  groujid  the 
crew  for  a  long  period  or  for  life.  Since  a  MOCOP  mission  could  involve  very 
senior,  highly  trained  people,  such  a  situation  might  be  very  undesirable. 

Taking  all  the  above  factors  into  account,  the  proposed  radiation  criteria 
is  based  on  the  philosophy  that  if  the  postulated  radiation  environment  docs  occur, 
the  crew  will  receive  a  dose  which  does  not  interfere  with  the  performance  of 
their  duties  but  which  may  result  in  some  minor  long-term  biological  effects. 

HIGH-ALTITUDE  ORBIT  TOLERANCE  CRITERION 

The  most  common  method  of  specifying  radiobiological  tolerance  criteria 
is  in  terms  of  cumulative  dose  (i.  e.  ,  the  arithmetic  sum  of  all  radiation  doses 
received  during  a  mission  should  not  exceed  a  given  amount).  For  short-duration 
missions,  this  is  certainly  a  valid  approach.  However,  for  extended  duration 
missions,  such  as  the  high-altitude  mission  under  consideration  here,  this 
method  may  not  be  realistic  and  may  lead  to  overly  pessimistic  shielding 
requirements.  For  long-duration  missions,  it  may  be  more  realistic  to  consider 
the  body's  ability  to  recover  from  sustained  radiation  damage.  A  criterion  formu¬ 
lated  on  this  basis  appears  to  lead  to  rehixed  shielding  requirement — either  in 
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terms  of  lower  shielding  weight  for  a  fixed  exposure  time  or  extended  exposure 
time  for  a  fixed  shielding  weight.  Application  of  this  criterion  to  the  high-altitude 
mission  is  discussed  below. 

Analytic  Formulation 

There  are  experimental  indications  that  biological  damage  is  influenced 

(15) 

by  radiation  rate  and  interval  effects. '  ^  Rate  effect  signifies  that  N  rad  of 

radiations  received  over  a  year  may  not  be  as  damaging  as  N  rad  received  over 

an  hour.  Interval  effect  signifies  that  two  applications  of  N  rad  spaced  months 

apart  may  not  be  as  damaging  as  two  applications  of  N  rad  spaced  several  days 

apart.  These  two  phenomena  seem  to  imply  that  there  is  an  ability  in  biological 

systems  to  repair,  at  least  in  part,  sustained  radiation  injury.  This  hypothesis 

flSl  ri7l 

is  formulated  analytically  by  Blair  ^  and  modified  by  Davidsom  ^  as 

r  1 

(t)  =  [f  +  (1-f)  e  ""  J  (13) 

where 

D  =  effective  dose 

e 

=  applied  dose 

f  =  nonrecoverable  fraction  of  injury 

t^  =  characteristic  recovery  time 

f  18l 

The  value  of  f  proposed  by  Blair  is  0. 10.  Schaefer^  ^  suggests  adding 

an  additional  0. 12  to  Blair's  value  to  account  for  the  effect  of  low-energy  solar 

protons  which  allow  no  biological  recovery.  The  value  of  t^  may  be  taken  as 

36.  1  days  in  accordance  with  the  recovery  half-time  value  for  man  of  25  days 

[17] 

suggested  by  Davidson. 
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K(|u:itioii  (13)  is  applicable  for  finding  the  residual  dose  at  time  t  due  to 
an  essentially  instantaneous  dose  applied  at  t  =  0,  and  hence  may  be  used  for 
calculating  the  dose  resulting  from  a  solar  flare  at  times  much  longer  than  the 
flare’s  duration,  or  the  dose  resulting  from  rapid  traversal  of  the  Van  Allen 
Belt.  Kquation  (13)  can  be  adapted  to  obtain  the  dose  due  to  a  prolonged  radia¬ 
tion  application  (such  as  cosmic  rays).  If  D  is  the  variable  dose  rate,  then 


t  -(t-T)/t 

f  +  (1-f)  e 


] 


D  (t)  dr 


(14) 


For  a  constant  rate  of  application  D,  (such  as  will  be  assumed  for  cosmic  rays), 

K 

i;(|.  (14)  leads  to 


D  (t)  =  D,  tf  +  t  (1-f) 
e  k  c 


Quantitative  Criteria 


(15) 


Although  there  is  a  range  of  reaction  in  people  subjected  to  the  same  acute 

radiation  dose,  it  is  generally  held  (e.g.  ,  Reference  19)  that  doses  in  the  range 

100  to  200  rem  produce  little  or  no  acute  radiation  sickness,  although  chronic 

[20] 

ellects  are  known  to  occur.  Also,  data  given  by  Pickering  (p.  42  of  Refer¬ 
ence  15)  indicate  that  doses  in  this  range  will  have  little  or  no  effect  on  an  air¬ 
crew's  efficienc}'. 

[2l| 

Baum  utilizes  the  recovery  concept  formulated  by  Blair  and  Davidson, 
and  described  in  the  previous  section,  in  proposing  a  radiobiological  tolerance 
criterion  for  astronauts.  He  suggests  that  it  is  permissible  for  an  astronaut  to 
receive  four  or  five  100  rad  doses  if  there  is  a  120-day  recovery  period  between 
the  exposures.  This  may  be  translated  as  a  total  of  375  to  400  rad  per  year  of 
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applied  radiation.  By  the  time  the  astronaut  is  ready  to  receive  the  fifth  100- 
rad  dose,  he  will  have  an  accumulated  residual  dose  of  about  40  rad  (see 
Eq.  (13)).  Hence,  the  effective  dose  at  the  fifth  exposure  is  100  +  40  =  140 
rad.  It  should  be  noted  that  Baum  expresses  his  criterion  in  terms  of  the 
absorbed  dose  unit,  the  rad,  while  we  have  been  using  the  effective  biological 
dose  unit,  the  rem.  Since  rem  =  rad  x  RBE,  where  RBE  ^  1  for  radiations 
considered  herein,  we  will  be  conservative  by  using  the  numerical  values  of 
Baum  in  terms  of  rem. 

The  present  study  modifies  what  is  essentially  Baum's  criterion  (with 
slightly  more  conservative  numerical  values),  and  takes  into  account  the  slowly 
applied  dose  resulting  from  cosmic  rays.  The  proposed  radiation  criterion 
for  the  long-duration,  high-altitude  orbit  assumes  that  the  effective  dose  compon¬ 
ents  are  defined  by  Eqs.  (13)  and  (15),  and  requires  that  the  absorbed  dose 
meets  all  three  of  the  following  conditions. 

(a)  The  instantaneous  applied  dose  shall  not  exceed  100  rem. 

(b)  The  sum  of  the  instantaneous  applied  and  accumulated  residual  dose 
shall  not  exceed  130  rem  at  any  instant. 

(c)  The  total  applied  dose  during  a  year  shall  not  exceed  350  rem. 

It  may  be  noted  that  this  criterion  untilizes  values  for  acute  and  chronic 
doses  that  are  higher  than  those  usually  suggested.  Although  adherence  to 
this  criterion  may  result  in  some  chronic  radiobiological  effects,  it  is  felt  that 
this  hazard  is  only  commensurate  with  the  other  hazards  of  military  space  flight. 
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LOW-ALTITUDE  ORBIT  TOLERANCE  CRITERION 


The  tours  of  duty  that  are  assumed  for  this  orbit  are  in  the  range  of  one 
week  to  three  months.  It  is  assumed  that  negligible  biological  recovery  occurs 
in  these  time  periods,  so  that  a  cumulative  dose  criterion  is  applicable.  The 
cumulative  dose  limit  is  taken  as  100  rem. 
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SECTION  V 


SHIELDING  CALCULATIONS 


A  detailed  mathematical  description  of  the  radiation  transport  through  a 
shield  material  is  given  in  Appendix  1.  This  analysis  considers  primary  protons, 
secondary  protons  and  neutrons,  and  bremsstrahlung.  The  analysis  is  used  to 
carry  out  computations  on  an  IBM  7030  STRETCH  computer  to  determine  the 
biological  dose  behind  different  amounts  of  shielding  material  for  the  different 
components  of  the  radiation  environment.  For  the  high-altitude  orbit,  the 
spectra  of  the  radiation  environment  components  given  in  Section  III  are  used 
directly  as  input  to  this  program.  For  the  low-altitude  orbit,  the  modifying 
effects  of  the  geomagnetic  field  on  the  free-space  radiation,  as  discussed  in 
Appendix  II,  are  fed  in  as  an  intermediate  step  in  the  computations. 
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SECTION  VI 


RESULTS  AND  DISCUSSION 


GENERAL 

The  results  of  the  shielding  computations  are  shown  in  the  dose  curves  of 
Figs.  7  through  10  for  the  high-altitude  orbit  and  Fig.  11  for  the  low-altitude 
orbit.  The  values  for  shield  thickness  (actually  in  terms  of  shield  weight  per 
unit  area)  in  Figs.  7  through  13  are,  of  course,  the  thicknesses  of  the  inner, 
storm-cellar  shielding  material  that  are  added  to  the  basic  cylindrical  shell. 

To  facilitate  the  calculations,  the  outer  cylindrical  shell  was  rejilaced  by  a 
spherical  shell  of  the  same  diameter  and  thickness.  Since  the  outer  shell  does 
not  contribute  materially  to  the  charged  particle  shielding  effectiveness,  this 


Fig.  7.  Dose  vs  Shield  Thickness  for 
May  10,  1959  (Type  A) 

Fla  re  (High-Altitude  Orbit) 
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Fig.  8.  Dose  vs  Shield  Thickness 

for  November  12-15,  1960 
(Type  B)  Flore  Sequence 
(High-Altitude  Orbit) 
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Dose  VS  Shield  Thickness 
for  February  23,  1 956 
Flare  (High-Altitude  Orbit) 
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Fig.  10.  Dose  vs  Shield  Thickness 
for  Passage  Through  the 
Van  Allen  Belt  (High- 
Altitude  Orbit) 


Fig.  1 1  .  Primary  Proton  Dose  vs 

Shield  Thickness  for  three 
Major  Solar  Flares  (Low- 
Altitude  Polar  Orbit) 


POLYETHYLENE  SHIELD  THICKNESS, 
GM/CM2 


35 
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Fig.  12.  Secondary  Dose  vs  Shield 
Thickness  for  November 
12-15,  1960  Flare 
Sequence  (High-Altitude 
Orbit) 


NOTE 

ABSCISSA  REPRESENTS  AMOUNT 
OF  SHIELDING  MATERIAL  ADDED 
TO  BASIC  1.5  GM/CM* 
ALUMINUM  SHELL 
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THICKNESS  OF  POLYETHYLENE  ,  GM/CM  2 


Fig,  13.  Effective  Dose  vs  Time  (High-Altitude  Orbit) 


36 


approximation  does  not  significantly  affect  the  results  and,  further,  lies  on  the 

conservative  side.  The  aluminum-foamed  polyurethane-aluminum  sandwich  wall 

of  the  outer  shell  was  replaced  by  an  equivalent  solid  aluminum  wall  of  thickness 
2 

1. 5  gm/cm  ,  since  the  weight  of  the  foamed  polyurethane  was  a  very  small  part 
of  the  total  weight  of  the  wall.  The  biological  dose  is  based  on  a  unit  target  at 
the  center  of  the  storm  cellar.  In  carrying  out  the  calculations,  all  protons 
that  penetrate  the  shield  with  energies  less  than  4  Mev  are  neglected  since  they 
can  be  stopped,  for  instance,  by  45  cm  of  cabin  air  at  7  psia  or  by  a  thin  layer 
(.  015  cm)  of  body  tissue. 

The  main  consideration  in  the  choice  of  shielding  materials  was  the 
effectiveness  of  a  given  amount  of  material  in  attenuating  protons.  Among  the 
elements,  hydrogen  is  by  far  the  best  in  this  regard,  and  carbon  is  second. 
However,  the  use  of  hydrogen  was  not  deemed  practical  because  of  storage 
and  containment  problems,  and  the  large  volume  required  because  of  its  low 
density.  Consideration  was  then  given  to  compounds  containing  large  proportions 
of  hydrogen,  such  as  water  and  polyethylene.  The  proton  range-energy  curves 
for  these  two  materials  are  almost  identical;  polyethylene  is  slightly  more 
effective.  This  slight  edge  in  effectiveness,  plus  simpler  problems  in  contain¬ 
ment,  seem  to  indicate  a  superiority  of  polyethylene  over  water  as  a  shielding 
material.  Thus,  the  proton  shielding  materials  used  in  the  high-altitude  orbit 
calculations  were  carbon  and  polyethylene  (with  specific  gravity  of  1. 0),  although 
for  all  practical  purposes,  the  results  of  polyethylene  are  the  same  as  those 
for  water.  The  calculation  results  shown  in  Figs.  7  through  10  indicate  the 
marked  superiority  of  polyethylene  over  carbon  in  this  application  so  that  only 
polyethylene  was  used  in  the  low-altitude  orbit  calculations. 
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The  dose  calculations  for  the  high-altitude  orbit  considered  the  contribu¬ 
tions  made  to  the  total  dose  by  secondary  protons  and  neutrons.  As  can  be  seen 
from  Figs.  7  through  10  and  Fig.  12,  these  contributions  are  small  for  shield 
thicknesses  in  the  range  of  practical  interest,  and  therefore  secondary  radia¬ 
tions  were  not  considered  in  the  low-altitude  orbit  calculations. 

HIGH-ALTITUDE  ORBIT 

The  dose  curves  for  the  various  components  of  the  radiation  environment 

are  shown  in  Figs.  7  to  10.  Figure  10,  which  shows  the  dose  received  in  going 

through  the  Van  Allen  Belt,  also  illustrates  that  polyethylene  is  ineffective  in 

attenuating  the  secondary,  bremsstrahlung  radiation  caused  by  the  trapped 

/  2 

electrons.  Note  that  the  addition  of  1  gm/cm  of  lead  to  the  inside  of  a  polyethy¬ 
lene  shield  reduces  the  dose  bj'  about  a  factor  of  3.  However,  the  rapid  passage 
through  the  Van  Allen  Belt  for  the  particular  transfer  trajectory  considered  here 
results  in  only  small  amounts  of  bremsstrahlung,  so  that  the  lead  is  not 
necessary. 

Figures  7  through  10  show  that  the  dose  due  to  secondary  protons  and 
neutrons  is  generally  small,  and  can  be  neglected  with  respect  to  that  due  to 
the  primary  protons.  For  illustrative  purposes,  the  variation  of  secondary 
proton  and  neutron  doses  with  polyethylene  shield  thickness  for  the  November 
12-15,  1960  flare  sequence  is  shown  in  Fig.  12.  As  mentioned  previously,  the 
dose  due  to  the  electrons  was  found  to  be  small  compared  to  that  due  to  their 
bremsstrahlung  so  that  only  the  latter  is  shown  in  Fig.  10. 

Figure  9  shows  that  the  dose  from  the  February  23,  1956  flare  is  sub- 

[6,  22] 

stantially  larger  than  that  reported  elsewhere.  The  difference  in  results 

for  this  flare  is  largely  a  consequence  of  differences  in  the  assumed  time- 
variation  of  the  energy  spectrum.  We  have  assumed  that  the  flare  decays  in 
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[6,  22] 

time  less  rapidly  than  have  other  authors,  and  hence  our  computed  total 

doses  are  larger.  However,  our  computed  initial  dose  rates  do  agree  with 
those  obtained  in  References  6  and  22.  Since  the  time-variation  of  the  spectrum 
is  of  necessity  formulated  from  measurements  at  only  two  times,  there  are 
obviously  many  ways  that  one  can  postulate  the  spectrum's  time-variation.  In 
view  of  this  uncertainty,  the  dose  results  for  this  flare,  as  obtained  both  here 
and  elsewhere,  should  be  used  with  caution. 


The  radiation  environment  is  related  to  the  mission  profile  and  is  assumed 
to  be  encountered  in  the  following  manner.  First,  the  Van  Allen  radiations  are 
encountered  during  outward  traversal  of  the  belt.  Next,  the  solar  flares  are 
encountered  with  the  timewise  phasing  shown  in  Section  III  while  the  vehicle  is 
on-orbit.  Also,  while  the  vehicle  is  on-orbit,  the  cosmic  radiations  are 
encountered  at  a  constant  rate.  Finally,  the  Van  Allen  radiations  are  again 
encountered  during  inward  traversal  of  the  belt.  This  mission  radiation  profile 
may  be  used  in  conjunction  with  the  dose  curves  of  Figs  7  through  10,  the 
assumed  cosmic-ray  dose  rate,  and  the  postulated  radiobiological  tolerance 
criterion  to  determine  the  shielding  requirements. 


Assuming  that  the  radiation  recovery  criterion  is  valid,  Eqs.  (13)  and  (15) 
may  be  applied  to  determine  the  effective  dose  as  a  function  of  time  for  differ¬ 
ent  thicknesses  of  a  given  shield  material.  Then,  by  trial-and-error,  the 
minimum  thickness  that  satisfies  all  conditions  of  the  criterion  can  be  deter¬ 
mined.  The  above  procedure  was  carried  out  (using  f  =  0.10  and  t  =  36.1 

days)  and  it  was  found  that  the  minimum  thickness  of  polyethylene  for  a  12- 

2 

month  mission  is  8.  5  gm/cm  . 

The  corresponding  plot  of  effective  dose  versus  time  is  shown  in  Fig.  13. 
It  is  seen  that  for  this  case  the  condition  that  determines  the  shield  design  is 
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the  instantaneous  dose  from  the  type-B,  November  12-15,  1960  flare  sequence. 

It  ean  also  be  seen  that  the  sum  of  the  instantaneous  and  residual  dose  at  any 
time,  even  if  the  type-B  flare  oeeurs  in  the  12th  month,  is  less  than  130  rem. 
Also,  the  arithmetic  sum  of  the  applied  doses  is  well  below  350  rem  under  all 
conditions.  Since  the  shield  is  designed  by  the  instantaneous,  100-rem  dose, 

the  minimum  shield  thicknesses  for  the  6-  and  3-month  missions  are  also  8.  5 

2  1 1  ^1 

gm/em  polyethylene.  It  is  interesting  to  note  that  if  one  uses  Schaefer's^ 

suggestion  for  nonreeoverable  fraction  of  0.22,  rather  than  Blair's^^^^  proposed 

value  of  0.  10,  these  results  are  still  valid. 

If  the  model  solar-flare  environment  is  modified  to  include  a  February  23, 

1956  type  flare,  the  shielding  requirements  are  changed  dramatically  In  this 

instance,  the  shield  is  designed  by  the  instantaneous  dose  from  this  flare,  and  at 

shielding  thicknesses  necessary  to  reduce  the  dose  to  100  rem,  the  doscs  Iroiu 

all  other  sources  become  almost  negligible.  Figure  9  shows  that  the  necessary 

,  2 

shielding  under  these  conditions  is  approximately  40  gm/cm  of  polyethylene,  an 
increase  by  a  factor  of  4.  7  over  the  previous  condition.  However,  as  discussed 
previously,  there  are  uncertainties  associated  with  the  data  from  which  the  dose 
curve  for  this  flare  was  calculated,  and  the  results  should  be  viewed  with  caution. 

LOW-ALTITUDE  ORBIT 

Dose  curves  for  the  three  flares  under  consideration  are  shown  in  Fig.  11 
for  the  ease  of  the  low-altitude,  geomagnetic  polar  orbit  and  for  both  simple 
and  modified  geomagnetic  cutoffs  (see  Appendix  II).  Only  the  dose  due  to  jirimary 
protons  is  shown,  because  the  shielding  calculations  for  the  high-altitude  orbit 
showed  the  contribution  to  the  total  dose  due  to  the  secondarj’  radiations  to  be 
negligible.  Also,  calculations  were  carried  out  only  for  a  polyethjdene  shield 
as  the  high-altitude  orbit  ease  showed  this  material  to  be  much  superior  to 
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carbon  for  radiation  shielding.  It  should  be  emphasized  that  the  results  associ¬ 
ated  with  the  February  23,  1956  flare  should  be  used  with  caution  because  of 
the  uncertainties  discussed  above.  Figure  11  shows  that  use  of  simple  dipole 
theory  gives  underestimates  of  the  dose  received  in  polar  orbits.  The  assumed 
modification  of  the  geomagnetic  field  due  to  the  magnetic  storms  accompanying 
the  solar  flares  is  found  to  increase  the  doses  by  23  to  30  percent. 

Comparison  of  the  solar  flare  dose  curves  of  Fig.  11  with  those  of  Figs. 

7,  8,  and  9  show  the  effect  of  the  geomagnetic  field  in  attenuating  the  proton 
dose.  This  attenuation  is  a  function  of  both  the  flare's  spectrum  and  the  shield 
thickness.  This  is  illustrated  in  Fig.  14,  where  the  ratio  of  geomagnetic-polar- 
orbit  dose  (using  simple  dipole  field)  to  free-space  dose  is  plotted  as  a  function 
of  shield  thickness  for  the  three  flares  of  interest.  These  results  show  that 
flares  with  a  "hard"  spectrum  (a  relative  abundance  of  high-energy  particles) 


Fig.  14,  Ratio  of  Polar  Orbit  to  Free-Space  Dose  vs  Shield  Thickness  for 
3  Major  Solar  Flares  (Using  Simple  Dipole  Geomagnetic  Field) 
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are  attenuated  less  bj'  the  geomagnetic  field  than  are  those  with  a  "soft"  spectrum. 
This  is  a  result  that  would  be  expected  since  relative  to  a  soft  spectrum  a  hard 
spectrum  has  a  higher  average  energy,  and  hence  the  effective  geomagnetic  cut¬ 
off  latitude  is  closer  to  the  equator.  This  results  in  the  vehicle  being  irradiated 
for  a  longer  time,  and  consequently  being  subjected  to  a  higher  dose  per  orbit. 
However,  it  should  be  noted  that  this  variation  of  dose  ratio  with  flare  spectrum 
hardness  is  not  significant,  and  for  a  given  shield  thickness  varies  by  only 
about  2  percent. 

The  components  of  the  model  environment  for  this  orbit  consist  of  one 

solar  flare  and  steady  cosmic  radiations.  Calculations  show  that  for  the  longest 

mission  under  consideration  for  this  type  of  orbit,  three  months,  the  cosmic- 

ray  dose  is  3  rem.  In  order  to  stay  within  the  cumulative  dose  limit  of  100  rem, 

the  solar  flare  must  result  in  a  dose  of  less  than  97  rem.  Assuming  that  a 

type-A  or  type-B  flare  have  an  equal  probability  of  occurring,  it  is  seen  from 

Fig.  11  that  the  shielding  requirements  are  dictated  by  the  type-B  flare,  and 

2 

that  a  polyethylene  shield  thickness  of  1.  0  gm/cm  is  required.  If  the  February 

/  2 

23,  1956  flare  is  considered,  the  requirement  jumps  to  5.  0  gm/cm  . 

TOTAL  SHIELD  WEIGHTS 

Up  to  this  point,  the  radiation  shielding  requirements  have  been  discussed 
in  terms  of  shielding  thickness  (in  weight/unit  area).  However,  the  quantity  of 
interest  to  the  mission  planner  and  preliminary  designer  is  the  total  weight  of 
radiation  shielding  that  must  be  carried.  Translation  of  the  shielding  require¬ 
ments  from  terms  of  thickness  to  terms  of  total  shielding  weight  is  straightfoi’ward 
if  one  stipulates  the  size  and  shape  of  the  space  that  is  to  be  shielded.  Basically, 
there  are  two  choices  available  to  the  designer:  to  shield  the  entire  vehicle;  or 
to  shield  a  minimum  size,  interior  compartment,  i.  e.  ,  a  "storm  cellar.  "  The 
former  alternative  is  probably  preferable  in  all  respects  except  a  very  important 
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one — weight.  For  example,  the  surface  area  of  the  vehicle  under  consideration 

2  2 
is  1170  ft  ;  the  weight  of  shielding  to  cover  the  entire  vehicle  with  8.  5  gm/cm 

of  polyethylene  (minimum  for  the  high-altitude  orbit)  is  20,  200  lb.  While  such  a 

weight  penalty  might  be  prohibitive  for  this  case,  the  weight  of  shielding  to  cover 

/  2 

the  entire  vehicle  with  1.  0  gm/cm  (minimum  for  the  low-altitude  orbit)  is  2390 
lb,  a  weight  penalty  which  might  be  acceptable. 

Two  factors  which  play  a  large  role  in  determining  which  shielding  concept 
will  be  used  are  the  total  weight  that  can  be  placed  in  orbit,  and  the  fraction  of 
this  orbital  weight  that  can  be  allocated  to  shielding.  The  former  depends  on 
the  capabilities  of  the  boosters  that  are  available  at  the  intended  time  of  launch¬ 
ing.  The  latter  depends  on  the  minimum  amount  of  equipment,  structure, 
propulsion,  etc.  ,  that  is  necessary  to  carry  out  the  mission. 

The  present  study  is  not  sufficiently  detailed  to  make  a  choice  between 

shielding  the  entire  vehicle  or  using  a  storm  cellar.  If  the  latter  course  is 

followed,  there  remains  the  problem  of  choosing  a  minimum  size  compartment, 

which  will  be  influenced  by  physical  and  psychological  factors,  as  well  as  by 

operational  factors  such  as  the  tasks  that  the  crew  must  carry  out  while  within 

the  shelter.  It  is  again  beyond  the  scope  of  the  present  study  to  propose  a 

storm-cellar  configuration,  but  if  the  crew  occupies  the  shelter  once  every 

few  months  for  one  or  two  days,  and  if  the  tasks  that  the  crew  must  perform 

while  within  are  minimal,  then  it  does  not  seem  unreasonable  that  a  spherical 

shelter  with  an  inside  diameter  of  7  ft  should  suffice  for  four  men.  The 

3  2 

enclosed  volume  is  then  179  ft  and  the  surface  area  is  154  ft  .  The  total 

2  2 

shield  weight  to  cover  this  shelter  with  8.  5  gm/cm  and  1. 0  gm/cm  are  2690 
lb  and  316  lb,  respectively. 
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It  should  be  noted  that  for  ease  of  shielding  calculations,  the  vehicle  shell 
was  represented  by  a  sphere  and  the  storm  cellar  was  assumed  to  be  an  inter¬ 
nally  concentric  sphere.  Also,  the  shielding  effects  of  any  internal  equipment 

[22,  23] 

were  not  considered.  It  has  been  shown  that  nonisotropic  distribution  of 

shielding  material  and  of  equipment  can  have  a  very  large  effect  on  the  internal 
dose.  However,  these  are  refinements  which  cannot  be  taken  into  account  in 
this  preliminary  study. 

Finally,  the  following  are  a  few  suggestions  that  may  be  considered  for 

reducing  the  shielding  weight.  One  possibility  for  reducing  the  total  weight  is 

through  the  use  of  nonhomogeneous  thickness  shields.  These  shields  can  be 

designed  to  provide  greater  protection  to  the  more  radiation-sensitive  parts  of 

the  body  such  as  the  head  and  torso,  and  less  protection  to  less  radiation- 

sensitive  parts  such  as  the  arms  and  legs.  Another  suggestion  is  to  use  the 

food,  water,  and  waste  materials  for  shielding.  Since  the  consumable  stores 

consist  mainly  of  hydrogenous  materials  of  high  shielding  efficiency,  in  amounts 

[24] 

sufficient  for  each  man  to  consume  10  Ib/day,  ^  they  should  be  an  effective 
adjunct  to  the  basic  shield  material.  The  internal  equipment  could  also  be  placed 
to  give  additional  protection  to  radiation-sensitive  parts  of  the  body.  It  is 
suggested  that  complete  reliance  should  not  be  placed  on  consumables  and  internal 
equipment  for  radiation  protection,  as  studies  have  shown  that  it  is  difficult  in 
practice  to  achieve  uniform  distribution,  and  that  "holes”  in  the  shielding  can 

[22  23] 

have  effects  out  of  proportion  to  their  size.  ’  In  times  of  acute  radiation 
danger,  such  as  an  encounter  with  a  flare  that  is  much  larger  than  expected,  a 
large  measure  of  protection  is  obtained  if  the  crew  huddles  together  and  thus 
use  their  group  self-shielding  ability.  Lastly,  for  large  space  stations  that  are 
proposed  for  the  somewhat  distant  future,  one  should  investigate  the  possibility 
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of  using  active  shielding  methods,  such  as  magnetic  fields  that  are  generated  by 
super-conducting  coils,  to  deflect  rather  than  stop  the  charged  particles. 

Dr.  Francis  W.  French 

fil .  KW"  P'  . 

Dr.  Kent  F.  Hansen 
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APPENDIX  I 


CALCULATION  OF  IU\DI/\TION  IX)SES 


PPLMARV  PROTON  DOSES 


The  integral  solar  flare  spectra  arc  characterized  bv  the  relation 


J  (>  R,  t,  ^)  =  (t)  exp 


R 


(Hi) 


where  -1(>R)  is  the  inteo'i’al  directional  flux  (protons/cm  -sec-stcr).  R  is 

Ihe  pi'oton  ripdity  (Mv),  P  is  a  unit  vector,  and  J  (t)  and  R  (t)  are  refer 

—  o  o 

enee  \  allies  foi-  the  integral  flux  and  rigidity,  The  reference  values  arc  func- 
lions  ol*  linio. 

The  pai'ticle  flux  per  unit  rigidity,  <I>  (R,  t,  P  ),  is  given  as 


<!>  (R,  t,  P) 


dj  (>  R,  t,  P  ) 

3r 


whereas  the  particle  flux  per  unit  energy  is 


<!•  (E,  1.  P)  -  ci>  (R,  t,  P) 


( 


protons 


cm  -scc-ster-Mev 


) 


A 


R  (t) 
o 


(E+E^,) 


^  2  ,,1/2 
[(E^K..)  -e;]  I 


(17) 
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The  proton  flux  is  assumed  to  be  uniform  in  the  space  between  the  earth's 
magnetic  field  and  the  sun.  Further,  the  flux  is  assumed  isotropic.  Consider 
the  spherical  shell,  an  octant  of  which  is  shown  in  Fig.  15. 


Z 


Fig.  15.  Shell  Geometry 


The  outer  radius  of  the  shell  is  a  (cm),  whereas  the  shell  thickness  is 
/  2 

A  (gm/cm  ). 

If  a  proton  with  direction  vector  f2'  is  incident  upon  the  shell,  then  the 
probability  that  the  proton  direction  vector  crosses  an  element  of  area  dA  at 
the  origin  is 


P  (S?  ,  n  )  dA  =  6  (SI'  +  n  ) 
—  — r  —  —  r 


2 

a 
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where  is  the  unit  vector  from  the  origin  to  the  outer  surface  of  the  shell. 
—  r 

The  probability  that  the  direction  vector  crossing  a  unit  area  at  the  origin  lies 
with  d  S2  of  n  is 


P  (n  ,  n)  df2  =  6  (n  +  n  )  dn 
' — r  —  —  — r  — 


The  directional  flux  of  protons  at  the  origin  of  energy  E  and  direction 

S2,  denoted  4>  (E,  t,  S2),  is  then  the  integral  over  the  spherical  surface,  over 

~  P 

all  directions,  and  over  all  energies  of  the  incident  flux, 

P 


(E,  t,  n)  =  y  dA^  dE'  y  dn*  $  (E',  t,  fi')  K  (E',  E,  A)^ 


P  n)  P  (£', 

where  K  (E',  E,  A)  is  the  probability  that  a  proton  of  energy  E'  penetrates 
a  thickness  A  and  emerges  with  energy  E.  The  surface  area  dA  is 


dA 


2 

a 


dn 

— r 


hence 


(E,  t,  n)  =  y  dE-  $  (E',  t,  2)  K  (E',  E,  A)  (18) 

The  function  K  (E',  E,  A)  is  given  as 

K  (E',  E,  A)  =  0  if  range  <  A 

K(E',  E,  A)  =  6[E-f(E',  A)]  if  range  >  A 
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r  o0l 

where  f  (E',  A)  is  found  from  the  range-energy  tables  for  protons.  The 
omnidirectional  flux  at  the  origin  is 

4>  (E,  t)  =  47r  4>  (E,  t,  S2) 

P  P  “ 

since  the  source  is  isotropic. 


The  dose-rate  density  in  ;  ~  at  the  origin  is  then 

gm-hr-Mev 


dD 

TPTTl 

— -  (E,  t)  =  Q  (E)  4>  (E,  t)  (19) 

at  p 


where  Q  (E)  is  given  as 


3  /sec\  “2  /rad\  -6 

Q(E)  =  3.6x10  j  10  jl.GxlO 


RBE  (E) 


=  5.  76  X  10 


RBE  (E)  g  (E) 


rem  1 


gm-hr  2 

proton/cm  -sec 


Thus  Q  (E)  is  the  dose  rate  per  unit  mass  per  unit  flux.  Values  of  the  RBE  (E) 
and  specific  energy  transfer  g  (E)  versus  E  are  given  in  Table  3. 

All  dose  rates  and  doses  reported  are  for  a  unit  volume  of  tissue  (H  0) 

2 

located  at  the  origin.  The  low-energy  limit  for  protons  was  taken  as  4  Mev, 
since  protons  of  lower  energy  would  be  stopped  by  the  air  with  the  capsule. 
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E  (Me 

4 

6 

8 

10 

15 

20 

30 

40 

60 

100 

140 

200 

250 

400 

600 

1000 

2000 

5000 

5100 


Table  3 


Dose  Rate  Per  Unit  Mass  Per  Unit  Flux 
(From  References  6  and  27) 


RBE  (E) 


rem 

rad 


g  (E) 


_ Mev _ 

gm/cm^-proton 


Q  (E) 


rem 

gm-hr-flux 


2.  6 

100 

1.  5  X  10  ^ 

2.  1 

72 

-2 

.  87  X  10 

1.7 

56 

. 55  X  lO”^ 

1.  5 

47 

. 41  X  lO”^ 

1.  3 

35 

.  26  X  10  ^ 

1.  2 

27 

.  19  X  10~^ 

1.  0 

20 

. 12  X  10~^ 

1.  0 

16 

.92  X  10~^ 

1.  0 

13 

.  75  X  lo”^ 

1.0 

7 

. 40  X  lO"^ 

1.0 

6 

.  35  X  10~^ 

1.  0 

5 

. 29  X  lO”^ 

1.0 

4 

. 23  X  10~^ 

1.  0 

3.  5 

. 20  X  lO"^ 

1.  0 

3.  0 

. 17  X  lO”^ 

1.0 

2.  5 

.  14  X  lO"^ 

1.  0 

2.6 

.  14  X  10  ^ 

1.  0 

3.  0 

.  17  X  lo"^ 

1.  0 

3.  0 

.  17  X  lO”^ 

51 


The  dose  rate  is 


dD 

rem 


(t) 


dt 


1 


dE  Q  (E)  4>  (E,  t) 

P 


while  the  dose  is 


D 

rem 


(20) 


(21) 


The  numerical  evaluation  of  the  dose  was  performed  as  follows:  a  discrete 

set  of  energies  E  was  selected;  the  flux  density  was  computed  by  use  of  the 
K 

range-energy  tables  and  Eq.  (17);  the  dose  rate  density  was  evaluated  from  Eq. 
(19);  and  the  dose  rate  was  computed  by  Eq.  (20),  using  the  trapezoidal  rule. 
Similarly,  the  dose  was  computed  by  Eq.  (21),  again  using  the  trapezoidal  rule. 

SECONDARY  PROTON  DOSES 

The  directional  flux  of  proton  secondaries  at  the  origin  is  denoted 


4.  (E,  t,  n) 


protons 

2 

cm  -sec-ster-Mev 


An  integral  expression  for  the  secondary  proton  flux  is 


4>^  (E,  t,  n)  =  ^dr^dE'^dn'  j  (£,  E',  D',  E,  f2)  j'dE"  Jd^' 

[  K  (E",  n";  E',  n')  2  (£,  E")  4>  (r,  E",  n")]( 

s  p  I 


(22) 


where  the  term  in  square  brackets  represents  the  birth  rate  density  of  second¬ 
aries  and  the  term  in  braces  represents  the  transmission  function  of  secondary 
particles. 
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In  particular 


T  (r,  E',  fl';  E,  fl) 
s 


probability  that  a  proton  at  position  r 
with  energy  and  direction  will 
cross  a  unit  area  at  the  origin  with  encrgj' 


E  and  direction 


H — 

\Mev-cm  -ster  / 


K  (E",  fi';  E',  fi')  =  probability  that  a  collision  by  a  proton 

S  ' 

of  energy  E"  and  direction  will 
produce  a  secondary  proton  with  energj' 

(protons  \ 

TT - i - I 

Mev-ster  / 

(£,  E")  =  macroscopic  proton  inelastic  collision 

cross-section  (em  ^). 


The  function  represents  both  cascade  and  evaporation  processes. 
The  evaporation  protcais  are  of  sufficiently  low  energy  (E'  <  25  Mev)  that  few 
penetrate  the  shield  and  may  be  ignored.  The  cascade  protons  arc  strongly 
directed  in  a  forward  cone  relative  to  the  primary  proton.  For  simplicity,  all 
cascade  protons  are  assumed  monodirectional  and  equal  to  the  primary  proton 
direction.  This  approximation  overestimates  the  secondary  dose.  However, 
ignoring  contributions  from  secondaries  produced  by  primaries  of  different 
directions  tends  to  compensate  for  the  approximation.  We  then  have 


(E",  n",  E',  fi')  =  n^  (E”,  E')  6  (fl"  -  fi') 


where  n  (E",  E')  is  the  number  spectrum  of  cascade  protons  produced  by 
s 

primaries  of  energy  E*^  The  cascade  proton  spectrum  is  simplified  by 
assuming  an  average  energy  for  the  secondaries.  An  empirical  relation 
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n  (E",  E')  =  n  (E”)  6 
s  s 


/e" -  -SI) 

\  . 445/ 


is  assumed.  Thus  we  have 


j^dE"  Jd^'  (E",  n";  E'  ,  fi')  2^  (r,  E")  4>  (r,  E",  fi") 

■  (rra)  (-  :®)  •  (*■ 


t 


protons 


m  -sec-ster-Mcv 


) 


The  transfer  probability  T  is  easily  derived  from  previous  results.  In 

s 


particular 


T  (r,  E',  f2";  E,  fi)  -  K  [E',  E,  A  (r)]  6  (fi'  +  )  —  a  (fi  +  fl) 

s  r  ^  r 

r 


Hence, 


4>^  (E,  t,  fi)  =  jr^  dr  jdfi^  jdE'  jdfi'  ^  a  (fi'  +  ) 

r"" 

*  a  (f2  +f2)n  2  (r,  -77^  W  (r,  — ^  ,  fi'N 

'  r  ^  s  y.445/  p  I  .445  /  y  .445  —  I 

=  jdE'  jdr  KlE',  E.  A  (r)l  ^  (r.  ■  «) 


The  omnidireetional  seeondary  proton  flux  at  the  origin  is 


(23) 


4>  (E,  t)  =  477  4>  (E,  t,  fi) 

s  s 
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The  numerical  integration  of  Eq.  (23)  is  performed  as  follows;  for  any 

energ>'  E,  the  various  values  of  E'  as  a  function  of  r  (i.e.  ,  A  (r))  are  found 

E' 

from  the  range-energy  tables;  the  corresponding  primary  energy,  _  ,  is 

.  445 


computed  and  used  to  find 


"s  (:^) ''p  f '  7^)  • 


the  range-energ>' 


tables  are  used  to  find  the  primary  proton  energy  outside  the  shell;  and  hence 
the  primary  flux  is  evaluated  from  Eq.  (17).  Having  found  the  secondary  flux 
at  the  origin,  the  dose  rate  and  dose  is  computed  as  before. 


The  data  used  for  the  calculations  is  given  in  the  following  tables. 


Table  4 


Number  and  Energy  of  Cascade  Proton 


Primary  Proton 
Energy  (Mev) 

Mean  No.  of 
Cascade  Protons 

Average  Cascade 
Proton  Energy  (Mev) 

50 

.56 

18 

100 

.  66 

36 

200 

.  86 

74 

400 

.  95 

148 

1000 

1.  12 

445 

2000 

1.  50 

890 

Table  5 


Proton  Inelastic  Cross-Sections 


[28] 


Proton  Energy  (Mev)  Inelastic  Cross-Section  (mb) 


25 

450 

50 

350 

100 

250 

200 

250 

400 

250 

1000 

250 
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SECONDARY  NEUTRON  DOSES 


The  directional  flux  of  neutron  secondaries  at  the  origin  is  denoted 


‘i>^  (E,  t,  fi) 


( 


neutrons 


cm  -sec-ster-Mev 


) 


The  integral  expression,  Eq.  (22),  holds  for  neutrons  with  the  appropriate 
modification  of  the  kernels  T  and  K. 

For  cascade  neutrons  we  again  assume  forward  scatter  and  a  primary¬ 
secondary  energy  correlation.  In  particular  we  assume 

(E”,  £?';  E'.  O')  =  a  (e"  -  a  (0"  -  O')  (24) 

where  the  subscript  nc  denotes  cascade  neutrons. 

The  transmission  function  for  neutrons  is  given  by 


T^  (r,  E',  n';  E,  fi) 


-T  (r. 


e 


E) 


6  (fi'  + 


2 

r 


6  +  n)  6  (E'  -  E)  (25) 


where  t  (r,  E)  is  the  optical  depth  of  the  material  for  neutrons.  Since  any 
collision  changes  the  neutron  direction  and  energy,  the  functions  e  '  ’  ^ 

and  6  (E  '  -  E)  appear.  Thus  we  are  computing  the  unscattered  neutron  flux 
at  the  origin.  The  contribution  of  scattered  neutrons  is  very  small  since  the 
mean  free  path  of  the  cascade  neutrons  is  considerably  larger  than  A. 
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The  flux  of  cascade  neutrons  at  the  origin  is  then 


The  numerical  integration  of  Eq.  (26)  is  performed  in  a  manner  analogous 
to  the  cascade  protons. 

Evaporation  neutrons  are  computationally  more  complex  since  the  evapora¬ 
tion  neutrons  arc  isotropically  emitted  from  the  nucleus.  The  birth-rate  density 
of  evaporation  neutrons  of  energy  E*  and  direction  is  then 


(27) 


The  transmission  function  for  the  evaporation  neutrons  is  the  same  as  for 
cascade  neutrons  and  given  by  Eq.  (25).  The  directional  flux  of  e\  apoi  ation 
neutrons  at  the  origin  is  then 


r  c"^ 

4>  (E,  t,  n)  -  \dr  - - - 

ne  —  J  Att 


(28) 


The  evaluation  of  the  angular  integral  in  Eq.  (28)  is  performed  by  a 
simplified  quadrature, 


(29) 
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where  /x  is  the  cosine  of  the  angle  between  the  proton  direction  and  the  radius 
vector.  Associated  with  any  value  of  p  there  is  a  fixed  shield  thickness.  The 
shield  thickness  and  range-energy  tables  are  used  in  conjuncticn  with  Eq.  (17) 
to  evaluate  the  directional  flux.  Equation  (29)  is  integrated  numerically  to  find 
the  omnidirectional  flux  of  protons  of  energy  E"  at  position  r. 

The  collision  cross-section  is  the  same  as  that  given  in  Table  5. 

The  function  n  (E",  E)  is  written  as 
ne 

n~~  (E",  E)  =  n~~  (E")  p  (E",  E) 
ne  ne 

with  n  (E")  the  mean  number  density  of  evaporation  neutrons,  which  is 
ne 

tabulated  in  Table  6,  and  p  (E",  E)  the  energy  distribution  function  for  evapora- 

f28l 

tion  neutrons.  The  evaporation  neutron  energy  spectrum  is  given  as^ 


p(E",  E) 


with  6  the  nuclear  temperature  which  is  correlated  with  the  incident  energy  as 


e 


2 


A 


(E") 


.  322 


with  A  the  atomic  weight. 

The  calculation  proceeds  in  a  manner  similar  to  the  cascade  neutrons. 
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Table  6 


Mean  Number  of  Cascade  and  Evaporation  Neutrons  per 

[  2  8l 

Inelastic  Collision 


Primary  Proton 

Mean  No.  of 

Mean  No.  of 

Energy  (Mev) 

Cascade  Neutrons 

Evaporation  Neutrons 

50 

,  30 

.  70 

100 

.  40 

.  70 

200 

.  60 

.  70 

400 

.80 

.  70 

1000 

1.  00 

.  80 

2000 

1.  20 

1.  00 

The  neutron  dose  in  rem  per  unit  neutron  flux  as  a  function  of  energy  is 
listed  in  Table  7. 


Table  7 


Neutron  Doses  per  Unit  Flux  as  a  Function  of  Neutron 

TT  [29] 

Energy 


Neutron  Energy  (Mev) 

12 

10 

8 

6 

4 

2 

1 

.  5 
,  1 


Dose  (Rem/hr-flux  x  10^) 

.  270 
.  270 
.  260 
.  230 
.  195 
.  148 
.  136 
.  083 
.  03 
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ELECTRON  DOSES 


The  Van  Allen  eleetrons  have  energies  ranging  up  to  ~  3  Mev  whieh  is 
sufficientlj'  low  that  all  eleetrons  will  be  stopped  in  the  outer  shell  of  the  capsule. 
The  bremsstrahlung  is  the  only  radiation  hazard  from  the  eleetrons. 

The  eleetron  omnidirectional  flux  is  denoted 


Only  1/2  of  the  flux  is  ineident  upon  the  spacecraft  surface  and  from  the  isotropy 
of  the  distribution,  the  eleetron  penetration  will  vaiy  from  0  to  the  mtiximum 
range.  Detailed  electron  trajectories  are  very  complieated  and  the  simplifying 
assumptions  are  made  that  all  eleetrons  stop  at  a  fixed  distanee  within  the 
shield,  and  further,  that  the  bremsstrahlung  is  produeed  at  the  point  of  stopping. 

The  angular  distribution  of  eleetrons  within  the  shield  will  approach 
isotropy  within  a  short  travel  distance  due  to  elastic  collisions.  It  is  therefore 
a  good  approximation  to  assume  the  bremsstrahlung  is  isotropically  distributed. 
The  omnidireetional  photon  flux  is  denoted 


Photons 


em  -see-Mev 


For  an  eleetron  of  energy  E,  stopping  in  a  material  of  atomic  number  z, 
the  bremsstrahlung  intensity  of  energy  hi^,  is  given  as^^^^ 


I  (hn)  =  2zk  (E-hi^)  per  electron 


with  I  (hn)  =  number  of  photons  times  the  photon  energ}- 
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Thus 


(hi^,  t)  =  4^  *  <i>  (E,  t)  0  <hi^  :<  E 

nP'  e 

E 

^  max 

=  ^  (E-hi^)  (E,  t)  dE  - |_2t2ilE -  (30) 

0  cm  -sec-Mev 

The  source  surface  is  located  as  in  Fig.  16. 


The  total  shell  thickness  is  A  .  The  thickness  of  material  from  the 

cm 

ETins 

source  surface  to  the  inner  edge  is  s  - — -  .  The  radius  to  the  source  surface 

cm 

is  r  (cm). 


The  probability  that  a  photon  emitted  from  the  source  surface  will  cross 

2  2 
a  unit  area  at  the  origin  is  1/4  tt  r  .  The  total  emitting  surface  area  is  47r  r  , 

hence  the  photon  flux  at  the  origin  is 
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^(hP,  t)  , 


'i>^  (ht^,  t) 


-T  (hi^,  s) 
e 


with  T(hi/,  t)  the  optical  path  for  photons  of  energ>’  ht^.  The  photon-dose- rate 
density  at  the  origin  is  the  intensity  at  the  origin  times  the  mass  absorption 
coefficient  for 


Mev 


dD(hp,  t)  -T  (hi^,  s)^  „  ,  a 

- ^ -  =  e  4>(hf^,  t)  hf^  — 

dt  p  see-gm-Mev 


The  dose  rate  density  in 


rem 


hr- Mev 


IS 


dD  (hi^,  t)  /u  X  ^  /u  X 

rem  ,  -T(hi^,  s)  ^  „  ,,  ,  a(hi^) 

- -  =  R(hi^)  e  4'  (hi^,  t)  hi^  — ' — 


=  II  (ht*)  e 


-T  (hj^,  s) 


(hi^) 


(hi^)  2kt 


hi^)  4>  (E,  t)  dE 


(31) 


The  dose  rate  and  dose  are  computed  by  integrating  Eq.  (31)  over  energy  and 
time. 

The  bremsstrahlung  dose  calculated  from  Eq.  (31)  represents  the  dose 
from  primary  photons,  i.  e.  ,  uneollided  photons.  The  build-up  of  secondary 
and  higher  order  photons  will  be  a  few  percent  at  best  and  is  ignored  in 
comparison  with  other  uncertainties. 
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APPENDIX  II 


POLAR  ORBIT  DOSES 


Below  the  earth's  geomagnetic  field  the  proton  spectrum  is  depleted  due 

to  deflection  by  the  geomagnetic  field.  The  cutoff  energy  is  a  function  of 

geomagnetic  latitude.  For  vertically  directed  protons  the  cutoff  rigidity,  P*, 
..  [31] 

is  approximated  as 


P*  =  P  cos  9 
o 


(32) 


with  6  the  geomagnetic  latitude  and  P^  =  14.9  where  c  is  the  speed  of 

light.  The  relativistic  relation  between  the  cutoff  rigidities  and  cutoff  cnergj’ 
E*  is  then 


E*  (9) 


2  8 
c  cos 


9  +  E 


(33) 


with  E^  (=  .  938  Bev)  the  proton  rest  energy. 


It  should  be  noted  that  Eq.  (32)  is  derived  by  assuming  that  the  geomagnetic 
field  is  a  simple  dipole  which  is  an  adequate  representation  in  times  of  low  solar 
activity.  However,  it  has  been  found  that  the  magnetic  storms  that  accompany 
solar  flares  act  to  decrease  the  cutoff  rigidity  at  a  given  latitude  from  that 
predicted  by  simple  dipole  theory.  This  results  in  an  effective  shrinking  of  the 
magnetically  shielded  region  around  the  earth  so  that  dose  calculations  made  on 
the  basis  of  simple  dipole  theory  give  results  that  are  larger  than  the  actual. 

An  average  correction  for  the  modified  geomagnetic  field  has  been  suggested  by 

[32] 

Webber  and  is  shown  below  in  Fig.  (17)  as  a  plot  of  effective  geomagnetic 
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latitude  versus  actual  geomagnetic  latitude,  for  use  with  Eq.  (32). 


50  55  60  65  70  75 

_ ^ _ J _ 1 _ I _ 1 _ I _ 

Effective  Geomagnetic  Latitude,  deg. 


45  50  55  60  65  70  75  80 

I _ I _ I _ I _ \ _ 1 _ I _ I 

Geomagnetic  Latitude,  deg. 


Fig.  17  Modification  of  Geomagnetic  Field  During  Solar  Flore 


[32] 


At  any  time  t,  a  vehicle  in  polar  orbit  will  be  shielded  by  the  earth's 
field.  Let  the  differential  proton  directional  flux  be  denoted  as 


<p{E  ,  a  ,  t) 


_ protons 

2 

cm  -sec-ster-Mev 


The  instantaneous  dose  rate  is  ,  at  time  t  , 

oo 

-  47r  r  dEK{E)  (^{E,  fi,  t)  (34) 

dt  J  '  '  '  _  sec 

E*(t) 

with  K  (E)  the  flux-to-dose  conversion  factor. 

The  integral  (Eq.  34)  may  be  divided  into  two  parts,  dose  due  to 
particles  at  or  below  the  maximum  cutoff  energy,  say  mid  dose  due  to 

particles  above  the  maximum  cutoff.  Hence, 
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m 


m 


m 


For  cosmic  rays,  the  proton  flux  in  free  space  is  assumed  constant  with 
time,  while  for  solar  flares  it  varies  over  the  life  of  the  flare.  Thus,  the 
computation  of  the  dose  due  to  solar  flares  involves  a  timewise  integration  over 
the  life  of  the  flare.  The  discussion  presented  herein  of  dose  calculation  in 
polar  orbits  is  in  terms  of  solar  flare  dose;  the  application  of  the  method  for 
calculation  of  cosmic  ray  dose  should  be  obvious. 

For  the  low-altitude  circular  orbits  being  considered,  the  vchick'^s 
angular  velocity  is  constant  and  equal  to  27r  radians  in  about  90  minutes.  It 
is  assumed  that  the  geomagnetic  field  is  symmetric  about  planes  containing  the 
geomagnetic  poles  and  the  geomagnetic  equator  so  that  for  a  constant  flux,  onl\‘ 
a  quarter  orbit  need  be  considered.  The  vehicle  traverses  a  quartej'  orbit  in 
about  22.5  minutes.  The  small  changes  in  the  solar  flare  flux  during  Uiis 
interval  are  neglected  and  the  flux  is  taken  as  constant  during  at  h'a.st  one 
quarter  orbit.  Therefore  an  average  dose  rate  during  a  quarter  orbit  is 
computed  as 


TT  /2 


m 


E*(0) 


o 


oo 


(3(1) 


m 
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=i  Itv; 


The  numerical  integration  of  Eq.  (36)  is  carried  out  as  follows:  A  discrete 
set  of  J  energies  E.  are  chosen.  For  E.  <  E*  the  corresponding  0,  (from 

j  J  ni 

j  J  o 

Eq.  33))  is  obtained.  Let  E.  =  E*  so  that  0.  =  0  .  For  energies  above 

J  m  1^ 

•'m  •'m 

cutoff,  the  contribution  to  the  dose  rate  is  obtained  as 


I  =  \  dE  K(E)  (p(E,fi,  t) 


AE^.  w.  K(E.)  (;P(E.,  n  ,  t)  (37) 


E* 

m 


J=J 


m 


with  the  weighting  coefficients. 

For  energies  below  cutoff  we  have,  (dropping  the  star  superscript) 


/ 

7r/2  m 

y  d0  y  dE  K(E)  <^(E,  n,  t) 
0  E(0) 


0 


90 


J  -  1  E 
m  m 

I  -I 

6,  E(0  ) 

■'m  •’m 


dE  K  (E)  <p(E,f2  .  t)  +  , 


0.  ,  E 
j-1  m 


0,  =  n/2  E 
1  m 


^  do  J  dE  K(E)(;p(E,f2,  t)  +  .  .  .  +y  d0  J  dE  K  (E)(;i?  (E,  fi,  t) 

0„  ~ 


0.  E(0.) 

J  J 


E(0,) 


(38) 


where 


m 


J  dE  K  (E)  cp(E,n,  t)  «  I 


E(..) 


m  “m 

i  <  '  "I 


(  ) 
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Define 


I  =  j  dE  K(E)  ip(E,  n.t)  (j  =  1,  2 . i^-j 


E(»j) 


Then 


(;!!D 


i  0=1  0^1 

E(0  )  E(0  )  E(8.^j) 


m 


(  )+...+ 


I 


(  ) 


E(0.  -  i; 
^ni 


=  I.  +  I.  +  . .  .  +  I. 

J  J+1  J^-1 

m 


(40) 


Therefore 


E 


m 

I 


dE  K(E)  cp(E,  n,  t)  «  - 


E(0.) 


(I.  +  I.  +  .  .  .  +  1.  -1)  +  (I.  -H  -1) 

J  J+1  jm  '  '  j-1  J  jjn 


'  J  W  ■ 


(41) 


Define 


-  1  -  \  I. 

S  =1  J 

j„-l  \  90  I  ; 
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(41') 


S. 

J 


-2 


m 


90 


+ 


S. 

J 


( 


0.-0 

J 


j±i 


90 


-) 


1 


) 


Then 


(43) 


of 


Calculation  of  the  total  dose  is  then  performed  by  numericjl  integration 


over  t  by  any  convenient  quadrature. 
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